Vortex lattice transitions in cyclic spinor condensates 
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We study the energetics of vortices and vortex lattices produced by rotation in the cyclic phase 
of F — 2 spinor Bose condensates. In addition to the familiar triangular lattice predicted by 
Tkachenko for 4 He, many more complex lattices appear in this system as a result of the spin degree 
of freedom. In particular, we predict a magnetic-field-driven transition from a triangular lattice to a 
honeycomb lattice. Other transitions and lattice geometries are driven at constant field by changes 
in the temperature-dependent ratio of charge and spin stiffnesses, including a transition through an 
aperiodic vortex structure. 
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One of the many remarkable properties of superfLuids is 
the appearance of vortex lattices in rotated systems [HQ. 
These lattices are periodic arrangements of vortices that 
allow the superflow outside the vortex cores to remain ir- 
rotational and are analogous to the mixed state of type-II 
superconductors in a magnetic field. Bose condensates of 
atoms with nonzero integer spin [H, HI, [H , referred to as 
"spinor Bose condensates" , combine spin and superfluid 
ordering in different ways depending on the spin and the 
interatomic interaction. These condensates, and the vor- 
tices and other topological defects that they allow, have 
been actively studied in recent years. 

Since the physics of individual defects in spinor Bose 
condensates is now understood for the most experimen- 
tally relevant cases with total spin F < 3, 0, 0,11, @, [ToL 
ll| . a natural next step is to understand physical situa- 
tions controlled by the collective physics of many defects. 
Two examples are the vortex lattice in a rotated conden- 
sate and the superfluid transition in a two-dimensional 
condensate. In general, the lowest-energy vortex defects 
of spinor condensates have both superfluid and spin char- 
acter. Although external rotation of the condensate cou- 
ples only to the superfluid part, the mixed nature of 
the vortices means that the interaction between the spin 
parts of different vortices is also important in determin- 
ing the vortex lattice. 

This letter uses a general approach to vortex lattice 
phases in spinor condensates, including the quadratic 
Zeeman anisotropy normally present in experimental sys- 
tems, to show that the cyclic phase of an F = 2 spinor 
condensate undergoes an unusual vortex lattice transi- 
tion in a weak applied magnetic field. This transition 
allows collective physics resulting from the nontrivial 
spin configuration of vortices to be imaged using spin- 
insensitive measurements. The comparison of energies 
of different lattices uses an Ewald summation trick that 
exactly reproduces previous results obtained for simpler 
lattices using elliptic functions More generally, the 
methods of this letter allow the energy of any periodic 



arrangement of vortices in a spin-anisotropic spinor con- 
densate to be rapidly calculated. We also show that un- 
der some conditions there is a strictly aperiodic vortex 
structure rather than a true lattice. 

Dilute F = 2 bosons interact via the potential V(\r\ — 
r 2 |) = S(ri -r 2 )(goPo+92P2+94P4,), where P F projects 
into the total-spin F state and gF = ^irTi 2 aF /M deter- 
mines <?f given of, the scattering length in the spin-i* 1 
channel. This two-body potential gives the interaction 
Hamiltonian [la, Il3j 
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with -0 a five-component vector field whose component 
ip m (r) destroys a boson at point r with F z = m, m = 
—2, . . . , +2, and F denoting the spin-2 representation of 
the SU(2) generators, ipt is the time-reversal conjugate 
of ip: iptm = {~^-) r ""4 } m- The parameters in this Hamilto- 
nian are determined by <7o,<?2, <74 via a = (3<?4 + 4^2 )/7, 
P = -(.92 - g4)/7, t = \{g - g 4 ) - j(g 2 - g±). To H int 
must be added the one-body Hamiltonian for an isotropic 
and spatially uniform condensate 



Ho 



dr y^V^t . VV> - ^H, (2) 



where fj,, the chemical potential. Minimizing this Hamil- 
tonian over single-particle condensates leads to three 
phases: ferromagnetic, antiferromagnetic, and cyclic. 
The cyclic phase that will be the focus of our work occurs 
when P, t > and is expected to be realized in a con- 
densate of 85 Rb atoms [l2j. The spinor structure of this 
state, having the symmetry of the tetrahedron, results in 
a nonabelian homotopy group which has been pursued in 
the liquid physics community for several years. 

In all existing experiments, an important effect even at 
the single-particle level is the existence of anisotropy in 
spin space resulting from the magnetic fields used as part 
of the trapping process. Including the hyperfine interac- 
tion, the bosons we consider interact with the external 
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FIG. 1: Orientations of the cyclic state in an external mag- 
netic field which breaks the spin rotational symmetry. Upon 
increasing the magnetic field, the spinor will undergo a tran- 
sition from state (a) to state (b). 



magnetic field as 7i z = VIS — 2[IbB z where T is the mag- 
nitude of the hyperfine interaction, I is the nuclear spin, 
S is the electronic spin, /is is the Bohr magneton, and B 
is the magnitude of the magnetic field taken to point in 
the z-direction. Within the manifold of spin-two states, a 
Hamiltonian which reproduces the correct energies up to 
a constant is given by Ti z = ^T 2 + (hbB) 2 + Y\ibBF z 
14j . This Hamiltonian can be expanded in powers of 



F z . Since the relaxation time of the total magnetization 
is typically longer than the condensate lifetime, the lin- 
ear term can be neglected. Particular attention has been 
paid to the next term which gives rise to the quadratic 
Zeeman effect Q. However, due to the high symmetry of 
the cyclic state, this quadratic term alone is not enough 
to select its orientation. For this case, one therefore must 
consider the cubic term which is at next order. 

To determine the spin states in the presence of such a 
magnetic field, one must also consider the spin exchange 
interaction energy of the condensate per particle which 
is E s = |n/3(F) • (F) + \nT\(x\Xt)\ 2 , where n is the con- 
densate density. Since the total spin is assumed to be 
conserved in experiments, we can neglect the first term in 
this expression. Minimizing E s + E z (where E z = (Hz}) 
over possible spinor states we find the following: At small 

magnetic fields, the spinor %ti = 0, 0, ^/2/3, 0^ 

(up to any rotation about the z-axis) is selected. In the 



classification scheme described in [15| , this state is repre- 
sented by a tetrahedron with one of its faces parallel to 
the xy plane. Upon increasing the magnetic field there 
is a transition at [IbE c = nr/16 to the spin orienta- 
tion xt2 = (sin(6l)/v / 2,0,cos(6'),0,-sin(6')/V2) T , where 
8 changes continuously with increasing magnetization. 
This spinor has the symmetries of a distorted tetrahe- 
dron with one of its edges parallel to the xy plane (when 
= -k/A it has the symmetries of the regular tetrahe- 
dron). These two types of orientations are summarized 
in Fig. [TJ The magnitude of the critical field is of the 



order of ambient fields in current experiments [l(| but 
smaller fields can in principle be simulated by optical 
means fl7| . 

Having identified the different types of tetrahedral 
states that are stabilized in an applied magnetic field, 
we now discuss the effects of rotation in addition to the 
applied field. The rotation couples to the phase of the 
condensate and has the effect of producing point vortices 
in two dimensions or lines of vortices in three dimensions. 
A vortex is a special type of configuration in an ordered 
phase that breaks a continuous symmetry: sufficiently far 
from a "core" region (linear in 3D or point-like in 2D) in 
which the order is destroyed. The configuration is locally 
in an ordered state, but cannot be smoothly deformed to 
the uniform configuration. 

We find in general that vortices form a two-dimensional 
lattice whose properties depend on the nature of the con- 
stituent vortices and the interactions between them. For 
simplicity, it is assumed that the magnetic field and the 
axis of rotation are in the same direction. Owing to the 
7r/3 spin rotation symmetry of the state (a) that is sta- 
bilized at fields B < B Cl its vortices are of three types: 
(n, m), (n — 1/3, to + 1/3) and (n + 1/3, to — 1/3), where 
n and to are integers and the first argument inside the 
parentheses is the winding number of the phase while 
the second is that of the spin. The vortex lattice that is 
formed has a net nonzero winding number for the phase 
and zero winding number for the spin. 

The energetics, in addition to the above constraints on 
the winding numbers, will depend on the stiffnesses K c 
and K s of the condensate corresponding respectively to 
the charge (phase) and the spin. (The expected behav- 
ior of these stiffnesses will be discussed in closing.) The 
interaction energy of two vortices (x\,yx) and (^2,2/2) 
separated by a distance r in the state (a) is given by 

E = 2nK c x 1 x 2 log(£/r) + nK syi y 2 log(£/r), (3) 

where £ is the typical radius of a vortex. For K s > K c , 
it is energetically favorable to produce only vortices of 
the type (1,0). However, for K s /K c < 1, the (1,0) vor- 
tex breaks up into (2/3,1/3) and (1/3,-1/3) vortices 
as can be seen by putting the appropriate values of the 
winding numbers into Eqn.[3] Once again, the energetics 
are subject to the constraints of rotation mentioned in 
the previous paragraph. For values of K s /K c < 1/4, 
each (2/3,1/3) vortex breaks up into (1/3,2/3) and 
(1/3, —1/3) vortices. Thus, in this regime, there are only 
(1/3, 2/3) and (1/3, —1/3) vortices with twice as many of 
the latter as the former. Since, the nature of the vortices 
is determined by the ratio K s /K c , so too is the lattice 
they form as will be described later. For B > B c , the 
state (b) is stabilized. This state can be shown to have 
vortices only of the type (n, m). When subjected to a 
rotation only (1, 0) vortices will bc produced like at low 
fields with K s > K c . 
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Thus, to summarize, the following kinds of vortices 
are produced by rotation: 1) (1,0) vortices for K s > 
K c or B > B c , 2) an equal number of (2/3,1/3) and 
(1/3,-1/3) vortices for 1/4 < K s /K c < 1 and B < B c 
and 3) twice as many (1/3,-1/3) vortices as (1/3,2/3) 
for < K s /K c < 1/4 and B < B c . In each case, the 
density of vortices is determined by the angular velocity 
of the rotation. 

Having determined the types of vortices that are pro- 
duced at the various values of magnetic field and stiff- 
nesses, we now evaluate the energies of vortex lattices. 
Due to the long-ranged nature of the logarithmic inter- 
actions, the energy of a vortex lattice is difficult to eval- 
uate directly. Thus, we develop a method that is simi- 
lar to the Ewald summation technique for the cohesive 
energy of three-dimensional ionic crystals 



18, m 



For 

simplicity, we use a scalar condensate to demonstrate 
the technique; the generalization to spinor condensates 
is straightforward and will be given presently. The en- 
ergy (in units of the stiffness) of a single vortex taken to 
be at the origin is given by 



R#0 \ / J \ 



(4) 



where R are the lattice vectors and £, the size of the 
vortex. The second term is due to a uniform nega- 
tive background charge of density po which arises from 
the fact we are working in a rotating frame of refer- 
ence. Note that each of these terms diverges individ- 
ually but their difference does not. The Ewald trick 
is to add and subtract a normalized Gaussian function 
-^pie~ r l<J from each point charge, where a is a screen- 
ing length. For instance, the potential of a point charge 
at the origin screened by such a Gaussian function is 
0(r) = ^Ei(r 2 /cr 2 ), where Ei(x) = die - * ft is the ex- 
ponential integral. Note that this will decay exponen- 
tially fast at large distances. In three dimensions the 
corresponding potential is <t>(r) — (1 — erf (r))/r where 
erf is the error function. 

Proceeding along these lines, the resulting potential 
corresponding to Eq. 2] is 
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(5) 



where G are the reciprocal lattice vectors and 7 is the 
Euler-Mascheroni constant. The first term comes from 
the density of point charges screened by the Gaussian 
function while second term comes from difference of the 
charge density of the Gaussian functions and the uniform 
charge density. The term — log ^Jt^ — J is obtained after 
subtracting off the additional Gaussian function at the 
origin where we omit the charge. Finally, the last term is 



to make the average of the screened potential zero [19(. 
On the other hand, the potential of a test charge away 
by d from the origin is 



0(d) = ^Ei 

R 



IR-d| 



E 

G^O 
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p0 2 a 
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The best check of this procedure is to see if the sum is 
independent of the parameter a. The two sums in real 
and reciprocal space in Eqns.[5]and[5]both converge ex- 
ponentially fast. In this way, the energy per vortex of the 

square lattice is found to be <f> s (0) = — log (j^ — 1.3105 
(where a is the lattice constant) while that of the trian- 
gular lattice at the same density is found to be </>t(0) = 
— log — 1.3211. Both are in precise agreement with 

the results obtained by integrating over the full spatial 
flow pattern [l[ . For vortices in spinor condensates, which 
contain windings of phase and spin, the above procedure 
is applied individually to each sector with the Ewald sums 
being weighted by the corresponding stiffnesses. The 
main advantage of the Ewald technique is that it can 
be generalized to treat complicated unit cells with an ar- 
bitrary number of vortices in them in a straightforward 
and numerically efficient manner. 

Let us first consider the case B > B c . As noted ear- 
lier, the vortices produced by the rotation are of the type 
(1,0). These form the usual triangular lattice for all 
values of K s /K c . For B < B Ci the fractional winding 
numbers of the fundamental vortices give rise to more 
interesting possibilities. For 1/4 < K s /K c < 1, the lat- 
tice is bipartite with equal numbers of (2/3,1/3) and 
(1/3, —1/3) vortices. We use the Ewald summation tech- 
nique to numerically evaluate the energy of the lattice 
assuming the same parallelogram unit cell for both sub- 
lattices and an arbitrary displacement between them. We 
then perform a minimization of the energy over these 
parameters to identify the lattice that is produced at 
different values of K s /K c . The sequence of lattices is 
described in Fig. [2j At exactly K s /K c = 1, the two 
sublattices do not interact with each other and each is a 
triangular lattice. As soon as K s /K c is lowered and the 
two begin to interact, the honeycomb lattice is stabilized 
and remains so till K s /K c = 0.76. Below this value, the 
vortices of one type move to the centers of the rhombic 
unit cells formed by the other type which we term an in- 
terpenetrating rhombic lattice. The internal angle of the 
rhombus changes continuously with K s / K c from 7r/3 at 
K s /K c = 0.76 to tt/2 at K s /K c = 0.64. The interpen- 
etrating square lattice thus obtained at K s /K c = 0.64 
is stable down to K s /K c = 1/4. This sequence of lat- 
tices is the same as obtained for rotating two-component 
condensates in the quantum Hall regime [jjcj . or equiva- 
lently the F = 1 polar condensate, but the values where 
the transitions occur are different for F = 2. 
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FIG. 2: A schematic depiction of the different types of vortex 
lattices obtained at different values of K s /K c and magnetic 
field. The color code for the vortices is: (1,0) yellow, (2/3,1/3) 
cyan, (1/3,-1/3) magenta and (1/3, -2/3) white. 



For K s /K c < 1/4, a lattice with (1/3,2/3) and 
(1/3, —1/3) vortices is obtained with twice as many of the 
latter as the former. Exactly at K s /K c = 1/4, the two 
sublattices do not interact and each is a triangular lat- 
tice. The sublattice of the (1/3, 2/3) vortices has a unit 
cell of length y/2 times that of the (1/3,-1/3) vortices. 
These two lattices are incommensurate for any rotation 
angle between them, which follows from showing that the 
nonzero squared lengths of lattice vectors in one lattice 
are disjoint from those in the other lattice. This incom- 
mensurability implies that the energy of interaction be- 
tween the two lattices can be calculated using the Ewald 
technique by averaging over all displacement vectors in- 
stead of specific lattice points, and the result is zero. In 
the other limit, K s /K c — ► 0, the interaction between all 
pairs of vortices is identical and a triangular lattice is ob- 
tained. While there are several way to distribute the two 
kinds of vortices in such a lattice, the lattice where the 
(1/3,-1/3) vortices form a honeycomb lattice while the 
(1/3,2/3) vortices are at the centers of each hexagon is 
the most symmetric one with three vortices per unit cell. 
The behavior between the incommensurate structure at 
K s /K c = 1/4 and this specific triangular structure as 
K s /K c — > is difficult to determine reliably by our tech- 
nique, since given the existence of the incommensurate 
structure, there is no justification for a numerical search 
over unit cells with a finite number of basis vectors. 

As demonstrated above, transitions between different 
vortex lattices can be tuned by a magnetic field B or the 
ratio K s /K c . While the field B can be applied directly or 
its effect simulated through optical techniques in experi- 
ments [l7j |. the ratio K s /K c is more difficult to manipu- 
late directly. In spinor condensates at low temperatures, 



this ratio is typically close to 1 but is renormalized by 
quantum and thermal fluctuations. Increasing temper- 
ature acts to reduce K s more rapidly than K c , because 
the soft spin modes that are excited at finite temper- 
ature have a larger phase space than the phase modes 
(assuming that the quadratic Zeeman term can be ne- 
glected). Both a nonlincar-sigma- model analysis and a 
study of Bogoliubov-like excitations suggest that under 
normal experimental conditions at nonzero temperature, 
K s is slightly less than K c so that the magnetic transi- 
tion will be observable. The best possibility to observe 
evolving vortex structure as K s /K c is further reduced is 
to raise the temperature very close to T c of the conden- 
sate: if the magnetic order is lost before the supcrfiuid 
order, as allowed by Landau theory, this ratio will rapidly 
decrease to zero in a narrow temperature range. 

To conclude, we have shown using the Ewald sum- 
mation technique that different types of vortex lattices 
can be produced in cyclic condensates as functions of 
magnetic field and the ratio of the charge and the spin 
stiffnesses. In particular, there is a magnetic-field-driven 
transition from a triangular to a honeycomb lattice. In 
the low-field limit, there are both abrupt transitions and 
continuous families of lattices as functions of the ratio 
of the stiffnesses, including the appearance of an incom- 
mensurate structure at one value. 
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